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We perform a systematic study of the odd-even mass staggering generated using a pairing 
interaction computed at first order in low-momentum interactions. Building on previous 
work including the (nuclear plus Coulomb) two-nucleon interaction only, we focus here 
on the first-order contribution from chiral three-nucleon forces. We observe a significant 
repulsive effect from the three-nucleon interaction. The combined contribution from two- 
and three-nucleon interactions accounts for approximately 70% of the experimental pair- 
ing gaps. This leaves room for higher-order contributions to the pairing kernel and the 
normal self-energy that need to be computed consistently. 

1. Introduction 

The nuclear energy density functional (EDF) approach is used to study medium- 
mass and heavy nuclei in a systematic manner Currently employed parameter- 
izations of the EDF provide a satisfactory description of low-energy properties of 
known nuclei. However, their empirical character and the spreading of the results 
obtained from different parameterizations, as one moves away from the valley of sta- 
bility and enters experimentally unexplored regions, point to the lack of predictive 
power of today's calculations. 
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Fig. 1. Different contributions to chiral 3N forces at N^LO. 

Our objective is to improve on this situation by developing non-empirical EDFs 
constrained explicitly from inter- nucleon interactions in free space. As a sta rting 
point, we performed the first systematic study of ground states of finite nuclei 
using a nuclear EDF whose^pairing part was obtained from low-momentum two- 
nucleon (NN) interactions ^^El first order. The present study builds on these 
results and focuses on the first-order contributions from chiral three-nucleon (3N) 
interactions. Details are reported elsewhere Higher-order contributions to the 
pairing kernel and the normal self-energy are left to future works. 

2. Method 

Our calculations are based on low- momentum interactions, consisting of an NN 
interaction Mowfc with a smooth cutoff A = 1.8 fm~^ and chiral 3N interactions Vs^ 
at next-to-next-to-leading order (N^LO) with Asn = 2.0 fm~ . The different 3N 
contributions are shown in Fig. [TJ For details and for the 3N couplings Ci^ cd and 
Ce used, we refer the reader to Ref. ^. 

Starting from the N^LO 3N force, we first construct an antisymmetrized, density- 
dependent two-body interaction Fnn by summing one interacting particle over oc- 
cupied states in the Fermi sea, extending the calculations of Ref. ^ to general isospin 
asymmetries. Fnn is then added to Mowfe to build the first-order pairing interac- 
tion from which the anomalous self-energy is computed. To proceed, we generate 
a high-precision representation of the pairing kernel as a sum of density-dependent 
separable terms and employ a local density approximation. Details will be given in 
Ref. 13. 

The remaining part of the nuclear EDF, i.e., the part accounting for the normal 
self-energy that describes the correlated single-particle motion, is taken as a semi- 
empirical Skyrme parameterization whose isoscalar and isovector effective masses 
have been constrained from Hartree-Fock calculations of symmetric and pure neu- 
tron matter employing the same NN and 3N low- momentum interactions The 
present results do not depend significantly on the Skyrme isoscalar effective mass, 
as long as the value of the latter at saturation density is 0.67 — 0.73. 

3. Results and conclusions 

We limit ourselves to discussing the odd-even mass staggering that provides a mea- 
sure of the lack of binding per particle in an odd isotope/isotone relative to its even 
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Fig. 2. Experimental and theoretical neutron/proton three-point mass differences along semi- 
magic isotopic/isotonic chains. Results are obtained from a first-order pairing interaction kernel 
with and without 3N contributions. 



neighbors; i.e., it essentially extracts the pairing gap Figure [2] displays theo- 
retical and experimental neutron/proton three-point mass differences along several 
semi-magic isotopic/isotonic chains from proton to neutron drip lines. Results ob- 
tained with and without 3N contributions to the first-order pairing interaction are 
compared. 

The main result obtained with NN only ^ is that theoretical neutron and proton 
pairing gaps computed at lowest order are close to experimental ones for a large set 
of semi-magic spherical nuclei, although experimental gaps are underestimated in 
the lightest systems. The addition of the first-order 3N contribution lowers pairing 
gaps systematically by about 30%. This is in line with the repulsive character of 
FsN in the ^Sq channel 1-^. Although the impact of the 3N contribution is rather 
insensitive to the structure of the particular nucleus under consideration, it displays 
a slight isovector trend along isotopic (isotonic) chains, most prominently seen in 
Fig. [2] for the neutron-rich lead isotopes 

The main conclusions of the present work are that (i) it is essential to include 
3N contributions to the pairing interaction for a quantitative description of nuclear 
pairing gaps, (ii) the first-order low-momentum results leave about 30% room for 
contribution from higher orders, e.g., from the coupling to (collective) density, spin 
and isospin fluctuations, (iii) in the next steps, the normal self-energy and higher- 
order contributions to the pairing kernel must be computed consistently starting 
from the same low- momentum NN and 3N interactions. 
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